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Minimal Flavour Violation hypothesis can provide an attractive framework for Dark Matter (DM). 
We consider scalar DM candidates carrying flavour quantum numbers and whose representation 
under the flavour group guarantees DM stability. They interact with the Standard Model fields 
through Higgs portal at renormalisable level and also to quarks through dimension-6 operators. We 
provide a systematic analysis of the viable parameter space for the DM fields, which are triplet 
of the flavour group, considering several DM-quark interactions. In this framework, we analyse in 
which cases the viable parameter space differs from Higgs portal models thanks to the underlying 
flavour structure. In contrast to minimal Higgs portal scenarios, we find that light DM in the GeV 
mass range as well as heavier candidates above Higgs resonance could be allowed by colliders, direct 
and indirect DM detection searches as well as flavour constraints. The large mass regime above mt 
could even be beyond the reach of future experiments such as Xenon IT. 



I. INTRODUCTION 

The determination of the nature of DM is one of the 
main challenges of cosmology and particle physics com- 
munity today. The Standard Model (SM) of particle 
physics, despite its great successes along with the likely 
recent unveiling of its last building block at LHC [1, 2], 
fails in providing a satisfactory DM candidate. In or- 
der to characterise the properties of the DM particle, we 
are forced to go beyond the SM. As of today, no clear 
scenario for New Physics (NP) is emerging. Concerning 
DM, we still ignore its intrinsic properties, how it couples 
to the SM particles and if the dark sector is made of one 
or several species. 

Within the SM, the SM scalar or Brout-Englert-Higgs 
("Higgs" for short) boson enjoys a special status since 
it can allow for a direct coupling to the dark sector at 
renormalisable level. It is already well known that the 
so-called "Higgs portal" provides a quite simple and at- 
tractive framework for DM phenomenology [3-20] . Inter- 
estingly, Higgs searches offer complementary bounds to 
direct and indirect DM detection searches on the viable 
parameter space of such DM models, especially in the 
low mass region [16, 21-28]. 

The DM stability requirement is usually ensured by im- 
posing by hand a discrete Z-i symmetry, under which the 
DM candidate is odd while the SM fields are even. Sev- 
eral works have however been investigating which models 
of NP could provide a more fundamental origin for its sta- 
bility, see e.g. Ref. [29] for a short review. A particular 
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remark should be make on models where more sophisti- 
cated global discrete symmetries [30-44] and horizontal 
gauge symmetries [45-48] have been adopted to deal with 
the DM phenomenology. 

Recently, Ref. [49] sketched the typical features of 
DM candidates coupling to quarks in a way that is 
consistent with the Minimal Flavour Violation (MFV) 
ansatz [50, 51]. Interestingly, the MFV, originally mo- 
tivated to suppress dangerous flavour changing neutral 
current (FCNC) processes in NP contexts, can have a 
dual purpose guaranteeing the stability of the DM. 

MFV is a working context that has been codified in 
Refs. [52-59] as a general framework built upon the 
flavour symmetry of the kinetic terms of the SM La- 
grangian. Focusing on the quark sector only, the latter 
presents a global flavour symmetry whose non-Abelian 
part is given by 

G f = SU(3) Ql x SU(3) Ur x SU(3) Dr . (1) 

Gf is only broken by the Yukawa interactions, unless 
the Yukawa matrices are promoted to be auxiliary fields, 
called spurions, transforming non-trivially under Gf. 
Without entering into details of a specific NP scenario, 
generic effects of flavour and CP violation can be de- 
scribed by means of an effective Lagrangian. If all the 
effective operators of this Lagrangian are constructed by 
means of the SM fields and the Yukawa spurions, all 
the flavour bounds are satisfied with a NP scale of few 
TeV [52-55, 58, 60-67]. 

A stable DM candidate can arise within the MFV con- 
text when the DM is assumed to carry flavour quan- 
tum numbers (see also Refs. [68-70]) and the DM field 
representation under Gf is chosen in order to prevent 
DM decay into SM particles. Following the approach 
of Ref. [49], we assume that the DM fields are neutral 
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under the SM gauge group, but transforms under Gf. 
For simplicity, we consider new scalar degrees of freedom 
as DM candidates with the lowest representation under 
Gf that guarantees their stability. By construction, the 
DM flavour-multiplets do not couple at the renormalis- 
able level to quarks. Their interactions with SM parti- 
cles result from the interplay between Higgs portal inter- 
actions and effective dimensions-six operators coupling 
DM to quarks. Notice that Ref. [49] proposed the rules 
for stable flavoured DM. In that paper, the DM analysis 
focused on scalar DM interacting through one particular 
d = 6 operator and with a fixed hierarchy in the DM 
components. The authors provided a survey of the typi- 
cal bounds that have to be taken into account from DM 
and flavour physics. This study missed though several 
new features in the DM viable parameter space that dis- 
tinguish the DM within the MFV context compared to 
the standard Higgs portal DM scenarios. 

In this letter, our purpose is to perform a systematic 
analysis of the viable parameter space of such DM con- 
text and to study how it differs from the simple singlet 
scalar DM interacting via Higgs portal. In particular, we 
obtain new features in the few GeV mass range, thanks 
to DM coannihilation induced by flavour symmetry, as 
well as larger viable parameter space for the DM mass 
range above the top mass niDM > wit, when an extra 
interaction to u— type quarks is considered. 

In the following, we first briefly summarise in Sec. II 
the MFV context and its extension involving a potential 
DM candidate. In Sec. Ill, we review the general form 
of effective d = 6 operators providing DM-quark interac- 
tions compatible with the flavour symmetry. In Sect. IV, 
we study the parameter space of this DM scenario con- 
sidering bounds from direct and indirect DM and collider 
searches and flavour constraints, and conclude in Sec. V. 



II. THE MFV-DM CONTEXT 

Deviations from the SM predictions induced by NP 
with generic flavour structure are already severely con- 
strained by the increasing accuracy in the determination 
of the CKM matrix elements, the measurement of a large 
number of FCNC processes and of CP asymmetries. In 
a general context, the scale of NP responsible for the 
flavour interactions should be above hundreds or thou- 
sands TeV [71] . On the other hand, under the MFV hy- 
pothesis which prescribes that all the sources of flavour 
and CP violation in any NP scenario are the same as 
in the SM, the scale of NP responsible for the flavour 
interactions is lowered down to few TeV. 

This is technically implemented by constructing effec- 
tive operators describing flavour and CP violation by 
means of quarks and Yukawa spurions, that transform 
under the flavour group Gf in Eq. (I). The SU(2)l dou- 
blet Ql and singlets u R and da transform according to 

Q L ~ (3,1,1), u R ~ (1,3,1), ^-(1,1,3), (2) 



while the Yukawa spurions transform as 

y„~(3,3,l), y d ~(3,l,3), (3) 

ensuring the invariance under Gf of the Yukawa La- 
grangian, 

& Y = -QLHY u u R -Q L HY d d R + h.c. (4) 

Quark masses and mixings are then correctly reproduced 
(but not predicted [57]) once these spurion fields get 
background values as 

Y u = y u , Y d = Vy d , (5) 

where y u ,d are diagonal matrices whose elements are the 
Yukawa eigenvalues, and V is the CKM matrix. 

Stable DM candidates can be found within the MFV 
framework looking for the representations of Gf that for- 
bid the construction of operators inducing the DM decay 
into SM degrees of freedom. It has been shown [49] that 
the lowest representation under Gf providing a stable 
DM candidate is a triplet under one of the SU(3)i com- 
posing Gf. The DM stability is then insured for any 
Lorentz representation of the DM candidate. For defi- 
niteness, we consider scalar DM S, neutral under the SM 
gauge group, and focus on the representation 

S~(3,l,l) (6) 

under Gf. Let us emphasise that our conclusions on the 
DM phenomenology should apply to any scalar triplet 
under one SU(3)i term of Gf, after a slight modification 
of the dimension-six operators that drive its interactions 
to the quarks (see Sec. III). 

The low-energy effective Lagrangian describing our 
setup is given by 

% = j?sm + d^s^d^s - v(s, h) + sejz + j?£ m , (7) 

where J^sm is the SM Lagrangian, V is the scalar po- 
tential involving the DM field S and the Higgs doublet 
H. Jzf e g could contain d = 6 pure-flavour operators, de- 
scribed in e.g. Ref. [52], that are suppressed by A^, while 
-5f c g M contains d = 6 DM-flavour operators suppressed 
by A 2 DM . Af and Ajjm are the characteristic mass scales 
of the messengers of the pure-flavour and DM-flavour 
interactions, respectively, and can a priori be distinct. 
While Af should be larger than a few TeV [52] to sat- 
isfy all flavour constraints, lower values of Ajjm are still 
allowed. In practice, we do not explicitly include the con- 
tribution of Jz?^ for the numerical analysis in Sec. IV C, 
as Jz? c L is not expected to affect the DM phenomenology 
given that in our analysis Ay is always larger than Adm- 
The SM scalar potential get modified at the renormal- 
isable level by the introduction of new scalar degrees of 
freedom. As already described in Ref. [49], in the MFV 
context, the new contributions to the scalar potential V 
read 

v d mis* (a i y + b (y„y„t) y + ...)S j+ 

+ A S* (a'Uj + b'(Y u Yj). l3 +...)Sj H^H, 
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where i,j are flavour indices, a, b, a', b' are dimensionless 
0(1) parameters and the ellipsis denotes further negligi- 
ble 1 Y d spurion insertions. After flavour and electroweak 
symmetry breaking (EWSB), (H) = v/y^2 with v = 246 
GeV, the S mass-squared matrix is diagonal and given 
by 



interactions. The coefficients c"j M take into account all 



J? m D~S* [m 2 A +m 2 B yl]S t 
where we have defined 



(9) 



f 



777 o<2 + — Av a' 



2 ' 



m% 



m 2 s b+^\v 2 b'. (fO) 



Of phenomenological importance are the trilinear cou- 
plings of the DM particles to the physical Higgs h, which 
in the mass eigenbasis read: 

&D~\vhSt (a' + b'ylJSi = ~\ivhS*Si, (11) 

where A l = A (a' + f/y« • )• 



III. THE DM-FLAVOUR OPERATORS 

S transforming as a triplet of Gf as in Eq. (6) pre- 
vents the construction of any operator containing a sin- 
gle DM field and quarks. A coupling between pairs of 
DM multiplets and quarks is still allowed through non- 
renormalisable interactions suppressed by A 2 DM . When 
the flavour symmetry Gf and the EW symmetry are un- 
broken, one can consider the following d = 6 operators: 



with 



cpDM 


1 5 

~~ A2 ^ c ijk£ {Og)ijM ' 
DM a=1 


(12) 






(13) 




= {u R ^u R] ){S*XSi) , 


(14) 




= (d Ri ^d Rj )(S* k %St) , 


(15) 


(°i)ijkl> 


= (Q Li u Rj ){S* k S e )H + h.c, 


(16) 




= (Q Li d Rj )(SZSt)H + la.c., 


(17) 



where i,j,k,£ are flavour indices (see also Ref. [72, 73]). 
The operators above correspond to an effective theory 
description of the NP sector giving rise to an additional 
sources of quark-DM coupling besides the Higgs portal 



possible flavour contractions and read as 
c \jki =cil«l« + c\lul kj + c 1 3 (Y u Y 1 l) lj l ke + 

+ c\l i:j {Y u Y^) H + cl{Y u Y^) u l kj + ... 
c %ki =c 2 lijl k e + cl(Y^Y u ) tJ l M + 

+ c \\ ld (Y u Yl) kl + c 2 4 (Yi)«(Y u ) kj + ... 

c%u =c?l«lw + 4{Y*Y d ) ij l ki + 

+ 4\ 13 {Y U Y^) M + cl(Yj) ie (Y d ) kj + ... 

c \jki =c l(Yu)ijlke + c 2^-ie(Y u )kj + 

+ 4(XMYuY^)kt + 4(Y u Y^) ie (Y u ) kj + 

c%u =4(Yd)ijlki + 4Ut(Y d ) k] + 

+ 4{Y d ) t0 {Y u Y^) u + cl(Y u Y^) u (Y d ) kj + 



(18) 
(19) 
(20) 

(21) 
(22) 



where we have considered all possible terms up to first 
powers in Y U Y^ and the dots refer to negligible contrac- 
tions associated to Y d insertions. In the following, we use 
real coefficients c", according to the MFV ansatz under 
which all the sources of flavour and CP violation should 
be associated to the Yukawas only. 

The scale Adm in Eq. (12) corresponds to a function of 
couplings and mediator masses whose exact combination 
depends on the UV completion of the effective theory 
under study. In principle, different operators or flavour 
contractions can depend on different type of messengers. 
01,2,3 would for instance be the result of a vector boson 
exchange, while ©4,5 would be due to the exchange of 
scalar or fermionic mediators. Here we assume that the 
structure of Eqs. (12)-(22) captures the DM phenomenol- 
ogy as long as the DM mass is lighter than the particles 
mediating the interactions. 

Also notice that when the energy scale involved in the 
physical processes is below the EWSB scale, both the 
gauge SU(2)l x U(1)y and the flavour Gf symmetry 
descriptions used in the definition of Eqs. (13)-(17) break 
down. One could then wonder for the appearance of new 
kind of operators, especially for those inducing DM decay, 
but this is not the case: the effective theory at low energy 
must match the theory at higher energies and therefore 
no new gauge or flavour couplings appear. 



IV. THE VIABLE PARAMETER SPACE 

In this section, we study the viable parameter space for 
DM analysing the impact of each dimension-6 operator 
Oi, ..,5 on DM phenomenology. 



1 Here we follow the common practice considering all interactions 
up to the first power in Y U Y^. See e.g. Ref. [61] for an anal- 
ysis studying the impact of resuming over the Yukawa coupling 
expansion in the context of MFV in the absence of DM. 



A. Constraints from DM and flavour physics 

Let us first summarise the constraints that have been 
imposed in each scenarios. We consider that A.; < ir 
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in order to preserve a perturbative regime. We take 
A; > — V^hj where Xh is the Higgs self coupling, to 
ensure that the scalar potential is bounded from below, 
assuming that the DM scalar self-couplings are up to 
~ 0(w). We also impose that Adm > mSi, m S2i m S3 
to avoid breaking down of the effective field theory de- 
scription. In addition, from DM and flavour physics, we 
have to take into account the following points: 

1. Relic abundance ^Idm- The total DM abundance 
deduced from cosmological observations constrains 
the thermally averaged effective annihilation cross- 
section (<tv) . In the simplest cases [74] , for the Tay- 
lor expanded annihilation cross section times the 



center of mass velocity av 



the Boltz- 



mann equations for the evolution of DM number 
density gives rise to the relation 



1.7 x 10 9 

y/g^(a + 3 6/: 



x f GcV- 1 
'xf)m p i 



(23) 



where h is the Hubble parameter 2 , Xf is the ratio 
between the DM mass and its freeze out temper- 
ature, g* is the number of relativistic degrees of 
freedom and m p i — 1.22 x 10 19 GeV is the Planck 
mass. The constraint floAih 2 ~ 0.1 [75] translates 
then into (av) ~ 3 x 1CP 26 cm 3 /s or equivalently 
(av) ~3x 10~ 9 GcV -2 . In our numerical analysis 
of Sec IV B, we use the code MicrOMEGAs [76, 77] 
that integrates more accurately the set of Boltz- 
mann equations and we impose 0.09 < floAih 2 < 
0.13. 

2. Direct and indirect detection: Direct detection 
searches are among the best test of Weakly In- 
teracting Massive Particle (WIMP) DM scenar- 
ios. This is especially the case of spin indepen- 
dent DM-proton scattering and for masses around 
50 GeV. The associated cross-section aoAip should 
be below the bounds of experiments such as PI- 
CASSO [78] and Xenon 100 [79] for m DM in the 
GeV-TeV range. Notice that such experiments 
have an energy threshold in the few GeV range and 
PICASSO's threshold is among the lowest ones. 
Indirect detection searches are also digging into 
the viable WIMP mass range. The annihilation 
cross-section times center of mass velocity av, can 
be tested by FERMI experiment [80]. The Cos- 
mic Microwave Background (CMB) observation 
experiments such as WMAP [75], SPT [81] and 
Planck [82] (in the very near future) provide com- 
plementary constraints on av. Indeed, the energy 
released into the Inter Galactic Medium (IGM) by 



2 Here, we keep the notation h for the Hubble parameter as it is 
widely used in the literature. It should not be confused with the 
physical Higgs field. 



DM annihilation can alter the thermal history of 
the Universe, leading to observable changes in CMB 
observables, see [83-86] and [87, 88] for the latest 
analysis. 

3. Colliders: Assuming that the particle resonance 
with a mass of about 125 GeV observed by the AT- 
LAS and CMS collaborations at the LHC [1, 2] cor- 
responds to the Higgs boson, one can constrain the 
invisible branching ratio to be Br(/i — > inv) < 0.15 
at 2a level [89]. Latest constraints from monojet 
events observed by the CMS [90] and ATLAS [91] 
collaborations can give rise to stringent bounds on 
a dm p complementary to direct detection searches 
for WIMP masses in the GeV range [92, 93]. Let us 
emphasise that the limits were derived for fermionic 
DM and in the case of scalar DM these bounds are 
not always that constraining. In the following, we 
use the results presented in Ref. [94]. 

4. Meson decays: The bounds on meson decays into 
invisible final state (see e.g. Ref. [73] for a re- 
view) strongly limit the direct couplings of quarks, 
apart from the top 3 , to the DM. In order to pass 
these constraints, one has to impose that JS-dm is 
larger than hundreds of GeV for low mass DM. 
Dimension-6 operators with such a large scale of 
NP can not guarantee the right DM relic abun- 
dance. Combining such bounds with the neces- 
sary small couplings A to avoid large Higgs invis- 
ible decay width, a DM particle coupling mainly 
to u — type and d— type quarks through 0\ 5 with 
mass rriDM < ™d/2 and ttidm < wi.b/2, respec- 
tively, is then excluded. 

5. Meson oscillations: Since no tree- level diagrams 
mediated by DM particles can be drawn, contri- 
butions to meson oscillations appear only at the 
loop level, as in the SM. The operators ©1,2,3 give 
negligible contributions to the meson oscillation ob- 
servables: the MFV ansatz ensures that the GIM 
mechanism also holds for these loop diagrams; an 
additional suppression oc 1 / A DM makes NP contri- 
butions smaller than the corresponding SM ones. 
In contrast, the operators O4 and O5 have a LR 
chiral structure and the corresponding contribu- 
tions to the meson oscillations are chiral enhanced. 
Following Ref. [49], the Wilson coefficients of the 
effective interactions (qR i qLj)(QR i Q.L j )-i for i > j, 



3 In Ref. [70] it was argued that single top production at LHC 
with large missing transverse energy involving flavour violating 
interactions can give rise to novel signature for DM detection. 
Notice though that for the DM models studied here, the t quarks 
only decay to u and c quarks and up to now LHC experiments 
can still not distinguish among light quarks jets. 
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depend on 



oc 



32^ 2 A 4 DM 



(V t *V tj ) 2 F 



' Q 

Oh 



(24) 



's, 



where ms^ms,) is the heaviest (lightest) of the 5 
components and F(x) — ((x + 1) loga;)/(a; — 1) — 2 
is the Inami-Lin function that vanishes in the limit 
x — > 1. For Adm > v, these NP contributions can 
be neglected while for Adm "C v, they can be large. 
Notice that for K meson system, the Wilson coef- 
ficients are suppressed by m 2 /A 2 DM and no sizable 
effects can be seen in AMr and e# (see Ref. [95] 
for a recent review for details on the meson oscil- 
lation observables) . Let us also mention that for 
Bd and B s systems, we expect new contributions 
in AMs d and AMb, only. Indeed, the absence 
of NP phases, in agreement with the MFV ansatz, 
prevents modifications on S^k s an d S^^. In addi- 
tion, let us stress that when the DM components 
are very near in mass the Wilson coefficient get an 
extra suppression through the Inami-Lin function. 
This is typically the case when the DM relic abun- 
dance is driven by coannihilations. 

For the models considered in this paper, only the 
DM-quark interactions through 04,5 can thus con- 
tribute substantially to meson oscillations. In prin- 
ciple, the operator O4 can give rise to non negli- 
gible contributions in D° — D° oscillations. The 
large theoretical uncertainties in the D meson sys- 
tem prevents though to set relevant bounds on C 4 
mediated interactions. This is not the case in B me- 
son systems. In the presence of O5 interactions, we 
impose that NP contributions from (24) are at most 
equal to the theoretical uncertainties on AMg d and 
AM Bs , i.e. AC ij = C% L /C^ M < 0.1. 

Let us also mention that, the operators O4 and C 5 
induce also modifications of the quark mass terms, 
when contracting the two DM legs at loop to draw a 
tadpole diagram: these contributions can be safely 
neglected, providing relative corrections suppressed 
by loop factors and by Tn DM /A DM . 



B. Analytical insights 

As previously discussed, the Lorentz structure of the 
operators 0\ ,..,5 can be categorised into two subgroups. 
The operators ©1,2.3 could arise from vector boson ex- 
change, while operators O4.5 would be associated to 
fermion or scalar exchange. At this point, it is instruc- 
tive to analyse the scale of NP Ajjm that would a pri- 
ori be necessary in order to give rise to the right relic 
abundance assuming a negligible contribution from Higgs 
portal. This is particularly relevant for DM masses be- 
low the Higgs resonance, where the bounds on the de- 
cay width of Higgs into invisible final state strongly con- 



strains the DM-Higgs couplings Xi (see the colliders sec- 
tion in Sec. IV A). 

In the limit of low center of mass velocity, v — > 0, 
the annihilation cross-sections associated to the processes 
S* Si — ► cjjqj for a fixed value of i and j become at leading 
order in v: 



N r c 2 v 2 m 



5m l m DM 



Am DM 



A&-kA du m 2 )M (1 - m 2 q /m 2 
N r c' 2 m 2 



.1/2 



DM J 



^ A DM 



(1 - m 2 q /m 2 DM ) 



3/2 



(25) 



where N c is the number of colors of the final state quark, 
c, d are a combination of 0(1) coefficients and moM and 
m q are the masses of Si and qj respectively. The first 
cross-section is p-wave suppressed while the second one 
corresponds to a s-wave driven process. Eq. (23) im- 
plies that for e.g. tudm = 50 GeV, Adm should be 
~ 450 GeV for 0\ (five families of quarks are involved) 
and ~ 200 GeV for O5 (mainly b quarks involved) in or- 
der to get the right relic abundance through annihilation 
only driven by the dimension-6 operators. Let us em- 
phasise though that given the velocity dependence of the 
<7i>'s above, it is clear that prospects for indirect detec- 
tion, involving velocities v ~ 10 -3 , will a priori be more 
constraining for ©4,5 than 01,2,3- bi practice, the coan- 
nihilations and Higgs portal interactions complicate the 
relic density analysis. The latter processes are fully taken 
into account in our numerical treatment of the DM mod- 
els with the MicrOMEGAs code after having introduced 
the proper Feynman rules using the LanHEP package [96] . 
We have also used MicrOMEGAs's tools for the calcula- 
tion of the cross-sections relevant for direct and indirect 
detection searches. 



C. Numerical analysis 

We have performed a random scan of the viable param- 
eter space for DM within the following parameter ranges: 



0.3 < |a, a', 6, &',<£| < 1-3 
0.1 GeV < m S3 < 1 TeV 
m DM < A DM < 10 TeV 

10~ 5 < |A| < 7T 



(26) 



considering real, positive and negative values for the pa- 
rameters a, a', b, 6', and A. Let us emphasise that once 
the value of ms a is fixed, one can obtain larger or smaller 
values of ms 1 and ms 2 by scanning over a,a',b,b' and 
A parameters. As a result, in contrast to Ref. [49], we 
have not imposed in our analysis that S3 is the lightest 
component of the S triplet. The constraints from per- 
turbativity, stability of the scalar potential, validity of 
the effective theory, and DM relic abundance are applied 
during the scan (see Sec. IV A for more details). 



6 




Figure 1: Operator 0\. Left: DM proton scattering cross-section as a function of the DM mass. Center: Annihilation cross- 
section as a function of the DM mass. Right: Scale of NP as a function of the DM mass. In all figures, the red points are 
excluded by the constraint on the Higgs invisible decay branching ratio, while blue (green) points (do not) pass the direct and 
indirect detection searches bounds. 



In Figs. 1 and 2, we present the results of the scans for 
the operators 0\ and ©4,5 in the planes mr>M — c^DMp, 
TfiDM — crv and m,DM — A-dm superimposing the con- 
straints from direct, indirect and colliders searches re- 
spectively. The corresponding plots for operators O2.3 
are very similar to the ones for 0±. The red points are 
excluded by constraints on the invisible decay branching 
ratio of 125 GeV Higgs, while the blue (green) points (do 
not) pass the constraints from the combined PICASSO, 
Xenon 100, FERMI and CMB experiments. 

Notice that in the itidm — av plots, we only show the 
constraint from FERMI on annihilation into bb at 100 % 
as they are the most limiting in the model under study. 
The limit associated to CMB represented in the plots re- 
sults from the analysis of Ref. [87]. The latter makes 
use of the recent SPT data assuming a value of the frac- 
tion of released energy in the IGM of 0.2 which should 
be representative of DM annihilating to quarks 4 . In the 
legend of the plots, the "viable" reference should be in- 
terpreted as "passes direct, indirect detection and CMB 
constraints". In the tudm ~ &DMp, we also present the 
reach of the future Xenon IT experiment [97] that will 
test most of the large mass regime of Or, ..,5 operators. 

We will see that some blue points do not pass mono- 
jet searches constraints from Tevatron or are within the 
reach of LHC at y/s = 14 TeV. In Figs. 1 and 2, we repre- 
sent with dotted lines the continuation of the limits from 
Ref. [94] to masses muu < 5 GeV. Notice however that 
in Ref. [94] slightly different assumptions were made on 
the coupling of DM to quarks compared to the present 



4 Notice that a similar analysis [88] shows the expected limits for 
each annihilation channels for WMAP7 and future Planck data 
release. It appears that the limits do not vary much going from 
annihilations into quarks to annihilation into gauge or Higgs bo- 
son. Also notice that Ref. [88] obtain a slightly less constraining 
bound for future Planck data. Here we stay a bit more conser- 
vative using the results of Ref. [87]. 



analysis. The lines in the mr>M ~ A-dm associated to 
Tevatron and LHC constraints are thus mainly a guide 
for the eye. We thrust though that Or, 2, 3 low mass can- 
didates are already excluded by such analysis (see Fig 1), 
while for O4 5 only a few models are within the reach of 
LHC. A more detailed analysis of LHC prospects for de- 
tection of such a DM candidate is however beyond the 
scope of the present paper. 



1, General features 

Several features appear in Figs. 1 and 2. Let us first 
concentrate on the av — m,DM plots. Typically, when the 
correct relic abundance is achieved and it is driven by an- 
nihilation processes, the annihilation cross-section today 
reaches ~ 3 x 10 -26 cm 3 /s. In this av range, we clearly 
see two preferred values of the annihilation cross-sections, 
especially in Fig. 2. The annihilation cross-section for 
Si. 2 DM today appears to be typically twice as large as 
the one for S3. 

This can be understood considering the fact that, in 
the MFV context, Si and S2 are generally quite degen- 
erate, see Sec. II. In the latter case, the effective DM an- 
nihilation cross-section [74] driving the relic abundance 
goes as (av) ~ (avu + avi 2 + av 2 \ + ci>22)/4 where av^ 
refers to the cross-sections times velocity associated to 
SiS* — > X annihilation and coannihilation processes (X 
is a SM final state) in the v — > limit, assuming that 
the relic abundance is driven by s-wave dominated pro- 
cesses. When coannihilation processes are suppressed, 
i.e. aviji^j <C avu and avu is identical for i = 1,2, we 
get (av) = avxi/2 ~3x 10~ 26 cm 3 /s. This is what we 
observe in e.g. the upper plot of Fig. 2 when Ajjm is large 
and more especially when the interactions driven by c 2 
terms in Eqs. (21)-(22) are suppressed. Notice that we 
apply the detection searches constraints as if all the DM 
in the galaxies was made of only one species Si- This 
is a rather conservative approach. If we had to consider 
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Figure 2: Operators O4 (left) and O5 (right). Top: DM proton scattering cross-section as a function of the DM mass. Center: 
Annihilation cross-section as a function of the DM mass. Bottom:Scale of new physics as a function of the DM mass. In all 
figures, the red points are excluded by the constraint on the Higgs decay invisible branching ratio, while blue (green) points 
(do not) pass the direct and indirect detection searches bounds. 



that half of the DM is made of e.g. Si and the other half of S2, some extra models with ms 1 < ms 3 could be 



considered as viable. 

Also notice that an important change in the number 
of relativistic degrees of freedom (from ~ 100 to ~ 10, 
see e.g. Ref. [98]) occurs in the early universe around 5 
Tqcd ~ 150 MeV due to the QCD phase transition. This 
implies larger values of the annihilation cross-section for 
rriDM ~ Tqcd x f . 

The relative absence of points in the m^M ~ mh/2 
is due to resonant DM annihilation through Higgs ex- 
change. This feature is also well visible in the Adm — 
rriDM and (TDMp — ttidm plots. Another dip appears 
in all cases at tudm ~ m\y when the annihilation into 
gauge bosons through Higgs exchange become kinemat- 
ically available. The other dips specific to the form of 
each operators are discussed in the following. 

Considering now the Adm — m-DM, we see that above 
the Higgs resonance the value of the Adm does not play 
an important role given that any value above 100 GeV 
is allowed by all the constraints. Below the Higgs res- 
onance, Adm should be typically smaller than (9(100) 
GeV to account for the DM relic abundance together 
with small contributions from the Higgs portal to evade 
a large correction to the Higgs width. Let us also remind 
that flavour observables, mainly meson decays into DM 
particles, impose that Adm is larger than hundreds of 
GeV for m DM in the GeV range [73]. More precisely, this 
constraint rules out tiidm < fnnjl an d rriDM < ttib/Z 
for the operators 04 and ©1,5, respectively and the ex- 
cluded regions are represented as (grey) shaded areas in 
Figs. 1 and 2. Moreover, the meson oscillation tends to 
constrain AC*-' < 0.1, see Sec. IV A. This threatens the 
DM models of O5 with the lowest values of Adm- The 
few points that are eventually excluded by this constraint 
are represented in light blue colors in the plots of the right 
column of Fig. 2. 

The candidates that pass direct and indirect detection 
searches constraints (dark blue points) appear in two dis- 
tinct mass ranges: below 10 GeV, where the direct de- 
tection bounds are weak, and above m/,/2 as in the case 
of Higgs portal models [16]: 

- For masses below 10 GeV, the scale of new physics 
is below ~ 200 GeV and the processes driving 
the relic abundance are typically related to the 
dimension-6 operators 0i,..,5- 

- Around 60 GeV, near the Higgs resonance, the scale 
of new physics is typically above 1 TeV and the 
coupling to the Higgs particle is small (A < 0.01). 

- Above 60 GeV, the scale of new physics is of several 
hundreds of GeV and the coupling to the Higgs 
particle is O(0.1 — n). 



5 A Tqcd ~ 150 MeV has been taken into account in MicrOMEGAS 
code 



2. Operators Oi,2,3 

In Fig. 1, we present scatter plots associated to the 
operator 0\. We have checked that equivalent results 
for 02,3 can be obtained (up to some 0(1) factor in the 
Adm)- Given the p-wave suppression of the annihila- 
tion cross-section, see Eq. (25), it is the direct detection 
searches sector that provides the strongest constraints 
when 0i drives the DM relic abundance. In the low mass 
regime, they are limited by PICASSO's energy threshold 
at vtldm ~ 4 GeV. Let us emphasise though that for 
scalar DM interacting through vector interactions with 
quarks, it has been shown [94] that monojet searches at 
Tevatron already provide very strong constraints on GeV 
range DM masses and rule out all the new potentially vi- 
able DM candidates with itldm < 10 GeV, as shown in 
the plot on the right hand side of Fig. 1. For masses just 
above Higgs resonance, the presence of extra interactions 
through 0i allows for lower value of the a dm p than in 
Higgs portal models and viable candidates that could be 
within the reach of future Xenon IT experiments and in 
some rare cases will be tested by LHC at y/s — 14 TeV. 

3. Operators 04,5 

In the case of the 4) 5 operators, the results are shown 
in Fig. 2. This time the annihilation cross-section is an 
s-wave process so that the indirect detection searches can 
compete with direct detection bounds for mDM < 10 — 20 
GeV. In the mDM — o~v plane, in addition to the Higgs 
and gauge bosons resonances, some other dips appear at 
mDM ~ m c ,m t in the O4 and mDM ~ TOb in the 05 
case. They are due to the fact that the main DM anni- 
hilation channels into S* Si — > cc, it and bb, respectively, 
are suppressed by factors of (— + vr? r)M ) x l' 1 in the low 
velocity limit, see Eq. (25). Away from those dips and 
resonances, the points getting lower values of the anni- 
hilation cross-section have their relic abundance mainly 
driven by coannihilation processes. 

As it can be seen in the bottom plots of Fig. 2, mono- 
jets searches [94] are not very limiting in the case of scalar 
type DM-quarks interactions. The low mass regime ap- 
pears to pass all constraints in contrast to the case of 
simple Higgs portal scenarios or to the case of Higgs por- 
tal + 0i.2,3 considered in the previous section. 

The particularities of operators O4 and 05 are the fol- 
lowing: 

a. Operator 04 : In the low mass range, viable DM 
candidates get an extra suppression of their relic abun- 
dance due to S1S2 — > uc coannihilation process, driven 
by the interaction terms proportional to c\ and c| in 
Eq. (21). For < m c , the coannihilation drives the 
relic density, while for 7711,2 > m c it is a combination of 
coannihilations and S2S2 — > cc that is relevant. Let us 
remind that in order to get a non-negligible contribution 
from coannihilations between two particles, their relative 
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mass difference should be Am/m ~ 0.1 in order to avoid 
exponential Boltzmann suppression. In the case of O4 
no particular tuning in mi 2 has to be invoked given that 
their degeneracy is granted thanks to the small u and c 
quarks Yukawa couplings in Eq. (9). 
Also notice that in the case of O4, for ttidm > TiH regime, 
one can obtain a very reduced value of <Jdm P compared 
to minimal Higgs portal scenarios (see e.g. [16]). This 
typically happens in the case of negligible coupling to 
the Higgs so that S. ( S* -!• WW, ZZ, it, hh through Higgs 
exchange are suppressed while O4 drives the DM relic 
abundance through SiS* — > it for ttidm > m t- This 
can be understood qualitatively by comparing the DM- 
proton scattering cross-sections for Higgs portal and O4 
models. In the first case cr£, Mp cx f 2 h X 2 /{m\m 2 DM )yL 2 [9] 

while in the latter case o- D \j p oc f^c 2 / (A\ )M m 2 DM )^ 2 
[49] and (J, is the DM- proton reduced mass. The nucleon 
form factors corresponds to f p h — 0.4 and f P 4 ~ 0.1 
in the Higgs portal and in the O4 case, respectively, us- 
ing the MicrOMEGAs default parameters [76]. The dif- 
ferences between f p h and f P 4 are principally due to the 
fact that O4 does not provide a coupling to the s quark. 
For large DM masses, one can also relates the values of 
av cx g 2 \ 2 /m 2 DM and av oc c 2 m 2 /Aj :iM for a DM abun- 
dance driven by Higgs portal or O4 processes, respec- 
tively. Eventually we get for tudm <C m t , a crjj\j p that 
is suppressed compared to <j\, Mp by a factor (f P 4/ f P h) 2 
and also by a factor ~ m\/{m 2 m 2 DM ). 



b. Operator 0$: Below 10 GeV, viable candidates 
have their relic density mainly driven by coannihilations 
Si S3 -> bd and S2S3 -> bs. All these coannihilation 
channels are described by the flavour interactions propor- 
tional to c\ and c| in Eq. (22). The S\$,S% coannihilation 
processes require some tuning of the parameters. Eq. (9) 
gives rise to m\,m<i ~ \m,A\ and 7713 ~ |ms| so that, if 
oCOjjCO = 77 l3 ~ 2m 12 . Allowing for different signs 
and the range of parameters in Eq. (26), we can get rela- 
tive mass differences ~ 0.1 necessary for coannihilation. 



From the above analysis appears that the coannihila- 
tion channels SiS^ ~* for O4 and Si S3* — > bd and 
S2S3 — > bs for C 5 enable viable DM parameter space 
in the GeV range. This is due to the flavour interac- 
tions driven by c% and cf with a = 4, 5 and this fea- 
ture is absent in standard Higgs portal models. No- 
tice that the GeV mass range will be tested in the next 
future by direct detection searches experiment, see e.g. 
Ref. [99] for SuperCDMS. Moreover, in the higher mass 
regime m dm > m t , O4 can give rise to rather suppressed 
DM-nucleon cross-sections that could even evade future 
Xenon IT constraints. 



V. CONCLUSIONS 



Embedding the DM problematic within MFV context 
can guarantee DM stability [49]. In addition, flavour 
DM-fermions interactions, that naturally arise in this 
framework, give rise to a richer DM phenomenology than 
in the minimal Higgs portal scenarios (see e.g. [16] for 
a recent analysis). In this paper, we focused on the 
MFV in the quark sector only with the Yukawas of the 
quarks being the only sources of flavour and CP vio- 
lation. For definiteness, we have restricted our analy- 
sis to the case of scalar DM fields, neutral under the 
SM gauge symmetry, but transforming as a triplet under 
one of the SU(2>) composing the flavour symmetry group 
G f = SU(3) QL xSU(3) UR xSU{3) dR . 

In this scenario, DM interacts with the SM fields 
through Higgs portal at renormalisable level and also to 
quarks through dimension-6 operators. We have consid- 
ered three vector operators and two scalar ones. We have 
performed a systematic analysis of the DM viable param- 
eter space for each of these operators. We have obtained 
that DM candidates can pass direct and indirect DM 
detection searches constraints both below 10 GeV and 
above m^/2. 

Complementary constraints from colliders and flavour 
physics have been taken into account. In general the non 
observation of meson decays into invisible final states ex- 
clude rriDM < "^zj/2 or hidm < mass ranges. 
In the case of dimension-6 operators with DM-quarks 
vector like interactions, monojet searches at colliders 
rule out all the new potentially viable DM candidates 
with moM < 10 GeV. For what concerns scalar type 
dimension-6 interactions, monojet searches are less lim- 
iting and allow for a viable rriDM < 10 GeV parame- 
ter space that will probably evade future LHC monojet 
searches. In most cases, it also evades constraints from 
flavour violating processes such as meson oscillations. In 
the large mass regime, most of the viable parameter space 
will be within the reach of future direct detection searches 
experiments such as Xenon IT, apart from the case of 
DM scalar-type coupling to u-type quarks that can evade 
the latter bound for mr>M > m t- 

To summarise, when the DM is embedded in the MFV 
framework, the DM stability is granted for certain rep- 
resentations of the DM fields under Gj. In this work, 
we systematically analysed the rich interplay between 
flavour and DM physics when considering vector and 
scalar type DM-quark interactions invariant under the 
original G/ symmetry. Compared with minimal Higgs 
portal models, DM candidates with ttidm ~ GeV range 
are viable in these scenarios. Furthermore, DM parti- 
cles with mrjM > m t could even be beyond the reach of 
future experiments such as Xenon IT. 
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